UNCLASSIFIED 

>0  419028 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIR6INIA 


UNCLASSIFIED 


ItOTICI:  Wien  government  or  other  drawings,  speci¬ 
fications  or  other  date  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related, 
foveraaent  procurement  operation,  the  U.  8. 
Qoveraaent  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever)  and.  the  fleet  that  the  Qovern- 
aent  nay  have  foraulated,  furnished,  or  in  any  my 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  Implication  or  other¬ 
wise  as  in  any  wanner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rljgits 
or  permission  to  aanufacture,  use  or  sell  any 
patented  invention  that  any  In  any  way  be  related 
thereto. 


Alt  UMIVEtSITY 
UNITED  STATES  Alt  FOtCE 


SCKHH.  W  ENGINEIR1NG 


MMUt-WHW*  M  WCI  MM, 


DO  C_ 

rpc'-^r-nr 

lV-  OCT8  '8S* 


a-WV^OCT  «  MN 


MBBUBPMPH 


Di  HI  JROON  RP  FUMA 
Lt  Ragraood  T.  Raurkat,  Jr. 
QM^Ftajr  a/63-3 


AF-wr-O-AUC  M  tr 


m  m  mam  bf  plasma 


TSBXS 

Presented  to  the  Faculty  of  the  School  of  Bnglneering  of 
the  Air  Force  Institute  of  Technology 
Air  University 

In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 


By 

Raymond  Thomas  Beurket,  Jr. ,  B.S. 
Lieutenant  USA 

Graduate  Nuclear  engineering 


August  1963 


(Ml/Pfays/63-3 


Th.-'  report  was  undertaken  mainly  for  the  purpoee  of  determining 
the  electron  energy  distribution  in  an  argon  17  excited  plasma  since 
until  the  writing  of  this  report,  no  energy  distribution  measure- 
ments  have  been  aade  in  any  17  plasma.  I  hare  written  this  report 
with  the  hope  that  a  person  with  a  general  physics  background  but 
with  little  or  no  knowledge  of  plamms  oan  appreciate  how  energy 
distributions  of  electrons  in  plasmas  nay  be  determined  and  specifi¬ 
cally  how  they  nay  be  node  in  17  excited  plasmas. 

Perhaps  the  most  valuable  aspect  of  this  investigation  (in 
addition  to  aoconpliahing  what  I  set  out  to  do)  was  the  disoovery 
that  no  natter  how  sisple  a  problem  sounds  on  paper,  unexpected  dif¬ 
ficulties  always  seen  to  crop  up  whan  now  ideas  are  tried.  There¬ 
fore,  la  scheduling  a  project  I  have  found  that  tins  Must  be  set 
aside  to  iron  out  any  unforeseen  difficulties  that  night  arise. 

In  conducting  this  investigation  I  received  valuable  assistance 
from  the  nanbers  of  the  llectronie  Technology  Laboratory  at  Vright- 
Patterson  Hr  force  Base.  I  an  indebted  to  Mr.  David  L.  Mays  for 
his  most  painstaking  explanation  of  equipment  and  helpful  suggestions 
and  to  Mr.  1.  Mason  Friar  for  the  use  of  his  radio  transmitter.  I 
an  also  indebted  to  Mr.  Varner  7.  Jahn  for  the  construction  of  the 
spherical  probe,  for  his  help  with  the  vaeune  qysten,  and  for  his 
humor  which  kept  up  the  author's  morale. 
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suggestions  he  made  throughout  the  experiment.  I  also  wish  to  express 
ny  appreciation  to  Professor  La  Verne  Lewis  who  was  my  thesis  adviser, 
and  to  Dr.  William  C.  Sppers  who  sponsored  the  project  and  made  the 
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ihrtnnti 

This  investigation  determined  that  electron  an ergy  distribution 
aeasureaents  in  an  XT  excited  pi  st  ara  faaaibia  If  tha  plii  is 
excited  by  a  cylindrical  coil  fad  push-pall  ft ran  a  balanced  power 
source  and  if  probe  neasursawnts  ara  aada  alone  the  axis  of  tha  coil. 
As  far  as  is  known,  tha  aaasareaants  aada  during  this  investigation 
are  tha  first  anargy  aaasaraasnta  aada  in  an  17  pi  a  was. 

Tha  baste  theory  of  plasaa  generation  and  general  probe  theory 
ware  reviewed;  Medicos'  spherical  probe  theory  was  used  in  asking 
the  determinations .  Kxperlaental  equipment  was  designed  and  built 
and  suggestions  for  laproveawnt  were  made. 

The  results  Indicated  that  in  the  argon  BF  plasaa  of  this 
investigation,  the  energy  distribution  is  markedly  non-Maxwallien. 
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BMumutnci  op  nuntoi  mmm  vmmmm 
n  ii  mujob  bp  FUAa 

i. 


Forty  years  ago  Lmagmlr  g ave  the  name  planaa  to  partly  ionised 
gaa  In  which  ttaa  positive  and  negative  oharges  ara  approadna tly  equal 
and  developed  the  baaio  plaana  theory  (Bef  8t3).  Sinoa  Langmlr'a  time, 
tba  field  of  plaaw  physics  baa  baan  expanded  greatly.  Planaa  theory 
it  aot  only  an  interesting  subject  but  la  iaportant  alnoa  probably 
99.9%  at  nattar  In  oar  universe  la  In  tba  pi  awn  atata  (Kaf  8i2) .  Thus, 
tba  study  of  planaa  la  lndaad  worthwhile. 

This  raport  daala  with  ona  araa  of  plaana  physios,  tba  datamlna- 
tion  of  alaotron  anargy  diatribation  In  an  BP  excited  plaana.  It  la 
important  baoauaa  tha  raanlta  of  this  experiment  nay  shad  light  In  other 
areas  of  investigation.  Tba  Electronic  Technology  Laboratory  at  Uright- 
Patteraoq  la  currently  doing  atudiaa  on  tha  gasaooa  laser,  la  an 
aactansion  of  this  thesis  topic  it  nay  be  possible  to  determine  tha 
electron  energy  distribution  in  the  laser  with  the  hope  of  using  the 
information  obtained  to  try  to  optimise  the  laser's  efficiency. 

It  present  no  work  has  been  done  in  making  energy  naaauranants  in 
an  BP  pi  asms,  and  tha  main  objective  of  this  raport  was  to  determine 
if  meaningful  naaauranants  could  be  made.  Measurements  ware  node  in  a 
cylindrical  argon  fillad  tube  with  the  restriction  that  the  measurements 
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war*  aad*  only  along  tfaa  cantor  11ns  of  th*  tub*.  Th*  reason  for  this 
raatriotion  will  bo  and*  apparent  In  Chapter  II. 


This  report  is  divided  Into  five  aaln  parts t 

(1 )  A  discussion  of  tho  fornation  and  proportlas 
of  a  pisses,  including  an  analysis  of  th*  con¬ 
ditions  under  which  aoasuraaants  can  be  aade 
in  an  17  plain; 

(2)  A  discussion  of  probo  theory  In  general; 

(3)  A  discussion  of  spherical  probe  theory; 

(4)  A  discussion  of  th*  oquipasnt  used  in  th* 

•xpor leant; 

(5)  Basalts  and  conclusions. 

A  ooaputor  pro  gran  which  was  used  in  reducing  th*  data  nay  b* 
found  in  Appendix  A.  Saai-log  plots  of  th*  second  derivative  of  the 
probo  curves  nay  bo  found  in  Appendix  B. 


2 


CBO/Pby  s/ 63-3 


ftfMfrlfiB  ABd  ffWKttM  SL  A  B— A 

When  «n  electric  field  is  applied  to  a  gas  at  ataospharie  pressure, 
the  gas  will  act  as  a  good  dieleotrie  unless  the  electric  field  is 
srtramsly  strong.  If  the  pressure  of  the  gas  is  reduced  sufficiently 
and  the  field  is  held  constant,  the  gas  will  become  conducting  and  the 
gas  wer  glow  with  a  color  characteristic  of  the  gas  employed  (Ref  1t132). 
In  this  experiment  argon  gas  was  used  and  the  resulting  plasaa  was  pink. 

The  field  applied  to  the  gas  any  he  either  a  DC  or  RF  field.  The 
DC  field  nay  he  generated  if  a  potential  is  applied  between  two  elec¬ 
trodes  placed  in  the  gas.  The  RF  field  nay  he  produced  by  surrounding 
the  gas  with  a  coil  and  then  applying  RF  energy  to  the  coll. 

Ionisation  hr  Rlectrons.  The  conducting  properties  of  the  gas  and 
the  glow  are  caused  by  partial  ionisation  of  the  gas  and  subsequent 
recombination .  This  partial  ionisation  of  the  gas  is  caused  by  the 
acceleration  of  electrons  by  the  electric  field  and  their  resultant 
energy  loss  by  inelastic  collision  with  gas  atoms.  The  primary  exciting 
electrons  nay  occur  initially  as  a  result  of  the  interaction  of  cosmic 
rays  with  the  gas  (Ref  14:27) .  If  this  naturally  occurring  process  is 
not  sufficient  to  produce  the  necessary  initial  source  of  electrons  for 
the  creation  of  the  plasma,  an  external  source  of  electrons  may  be  pro¬ 
vided  (as  was  the  case  in  this  experiment)  by  placing  a  Tesla  coil  near 
the  tube  containing  the  gas.  Once  these  primary  electrons  have  been 
provided  and  the  pressure  of  the  gas  and  applied  electric  field  strength 
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have  bam  proparlor  adjusted,  the  plaana  should  beccne  self-sustaining 
since  aaoh  ionisation  will  prodooe  at  laaat  ona  additional  alaotron 
eapabla  of  oaoaing  fnrthar  ioniaatlona  until  a  balance  between  genera- 
tion  of  alaotrona  and  losa  of  alaotrona  ia  achieved. 

It  ia  vail  known  that  a  partiela  of  oharga  q  whan  falling  throng 
a  potential  V  will  achieve  a  kinetic  energy  qV.  m  order  to  nore 
fully  undaratand  the  nethod  by  which  ionisation  take  a  place,  it  ia 
naoeaaary  to  undaratand  how  tba  alaotrona  loaa  their  kinetic  energy. 

Upon  achieving  a  certain  kinetic  energy,  the  alaotrona  will  proceed  to 
loaa  their  energy  in  four  principal  ways  (Ref  14s16)t 

(1)  The  alaotrona  surrender  energy  upon  suffering  collisions 
at  the  vail  of  the  tuba  and  to  the  anode  in  a  DC  dis¬ 
charge. 

(2)  The  electrons  transfer  energy  to  the  gas  in  the  fora 
of  heat  aa  a  result  of  elastic  collisions. 

(3)  The  electrons  excite  the  gas  atona  as  a  result  of 
inelastic  collisions. 

(4)  The  electrons  ionise  the  gas  atona  after  colliding 
inelastically. 

In  addition,  in  an  RP  plasna,  those  electrons  whose  direction  of  notion 
is  against  the  electric  field  will  also  lose  kinetic  energy. 

Since  gas  atons  of  atonic  weight  A  have  1840  A  the  ness  of  an 
electron  (an  argon  atoa  has  73,600  tljws  the  mass  of  an  electron),  the 
energy  loss  by  the  electrons  upon  elastic  collision  with  gas  atons  will 
be  so  snail  aa  to  be  negligible.  This  can  be  shown  to  be  true  if  the 
equations  of  conservation  of  energy  and  ncettntua  are  considered  and  the 
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electron  and  atom  are  assumed  to  be  —ell  eleatlo  spheres.  Thus  even 
though  the  electrons  suffer  elastic  collisions  In  falling  through  a 
potential  V  their  kinetic  energies  will  still  be  qV  to  a  high 
degree  of  aeeoraoj  (Ref 

If  the  energy  of  the  electrons  is  greater  than  the  energy  equiva¬ 
lent  of  the  ionisation  potential  of  the  gas,  the  electrons  nay  ionize 
the  gas  upon  colliding  inelastlcally  with  a  gas  atom.  If  the  energy  of 
the  electron  is  less  than  the  energy  equivalent  of  the  ionisation  poten¬ 
tial,  the  electron  will  lose  energy  upon  an  inelastic  collision  with  a 
gas  aton  and  will  excite  the  atom.  As  long  as  the  energy  gained  by 
the  electron  is  less  than  the  lowest  excitation  level  of  the  gas  atom, 
all  collisions  wist  he  elastio  and  if  all  the  electrons  have  energies 
that  low,  no  plasma  can  exist. 

The  electric  field  strength  and  pressure  mat  be  correct  to  insure 
a  plasm  because  the  potential  through  which  the  electrons  fall  per 
collision  is  the  factor  that  determines  whether  they  will  achieve 
enough  energy  to  cause  excitation  or  ionisation.  The  average  distance 
between  collisions  is  known  as  the  mean  free  path  of  an  electron  in  the 
gas,  and  its  symbol  is  X.  If  the  electric  field  strength  is  symbolized 
by  £,  then  the  potential  difference  per  mean  free  path  is  IX.  Since 
X  represents  the  average  distance  between  collisions,  if  the  number  of 
particles  increases,  then  X  must  decrease.  This  is  equivalent  to 
saying  that  as  the  pressure  rises,  the  mean  free  path  drops.  Thus,  there 
is  an  Inverse  relation  between  mean  free  path,  X,  and  pressure  pj  hence, 
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the  potential  difference  per  mean  free  path  will  be  proportional  to  E/p 
(Ref  14:16).  In  order  to  achieve  excitation  or  ionisation,  the  ratio 
Vp  *u*t  be  adjusted  so  that  it  is  large  enough  to  insure  that  the 
potential  per  mean  free  path  is  large  enough  that  the  electron  will 
gain  enough  energy  after  several  aean  free  paths  to  excite  or  ionise 
the  ate*.  This  wqr  be  accomplished  by  increasing  the  field  strength 
or  lowering  the  pressure.  In  this  experiment,  the  pressure  was  set 
at  0.2  torr  which  was  sufficiently  low  for  the  field  intensity 
involved  to  give  a  dense  plasma. 

Ianiartlffl  la  I&BA  ifiA  Photons,  it  this  point,  section  will  be 
made  of  the  positive  ions  and  their  roles  in  causing  ionisation. 

Positive  ions  are  formed  whenever  an  electron  is  ejected  from  a 
neutral  atom.  These  ions  will  also  gain  kinetic  energy  as  a  result 
of  being  accelerated  by  the  electric  field.  The  masses  of  the  ions 
are  essentially  the  same  as  the  masses  of  the  neutral  gas  atom  since 
the  loss  of  an  electron  will  affect  their  masses  negligibly.  Since 
the  ions  have  about  the  same  mass  as  the  neutral  gas  atoms,  the  ions 
will  lose  most  of  their  energy  after  each  collision  and  hence  in  most 
cases  never  gain  enough  energy  to  excite  or  ionize  the  gas  atoms.  Thus, 
we  can  essentially  ignore  the  effect  of  excitation  or  ionization  caused 
by  ions  (Ref  14:18). 

Electrons  may  also  be  liberated  by  the  photoelectric  effect,  but 
this  contribution  is  very  minor  too.  The  wavelength  of  the  photons 
must  be  short  enough  to  cause  excitation  or  ionization  and  this  is 
rarely  the  case.  The  principal  means  of  ionisation,  therefore,  is  by 
the  Inelastic  collision  of  electrons  with  gas  atoms. 
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Since  the  noaber  of  positively  end  negatively  charged  particles  in 
the  plasma  is  essentially  the  sane,  the  net  space  charge  in  the  plaana 
region  is  zero.  In  a  DC  plaana  the  potential  drop  of  the  tube  occurs 
across  the  cathode  and  the  positive  ion  sheath  that  foraa  very  near  it. 
The  rest  of  the  plasma  is  essentially  field  free.  This  is  fortunate 
because  when  a  plasma  is  probed,  the  probe  asst  be  in  a  field-free 
region  in  order  to  make  meaningful  measurements.  In  the  RF  case,  the 
electric  field  will  penetrate  the  plaana  and  the  field  lines  will  be 
concentric  about  the  axis.  Probe  measurements  are  not  impossible  in 
the  RF  plaana  if  the  probe  is  restricted  to  the  axis  of  the  tube  and 
a  push-pull  feed  is  used  to  excite  the  plaana.  With  this  arrangement 
there  will  be  no  RF  potential  between  probe  and  plaana  and  the  probe 
will  lie  in  a  field-free  region  as  the  following  discussion  indicates. 

If  the  plasma  is  generated  by  feeding  push-pull  a  coil  surrounding 
the  plasma  tube,  two  Important  conditions  occur  in  the  tube.  First,  a 
ground  plane  will  exist  in  the  tube  if  the  electrical  center  of  the 
coil  is  grounded.  Secondly,  at  the  geometrical  center  of  the  tube  the 
RF  electric  field  should  be  zero  if  the  exciting  frequency  is  not  too 
high. 

r.T»<mnri  PI  ma .  Feeding  the  coil  push-pull  implies  that  the  voltage 
impressed  on  one  side  of  the  coil  will  be  180°  out  of  phase  with  the 
voltage  impressed  on  the  other  side.  If  the  voltages  are  balanced  and 
if  Z  (Fig.  1,  shown  on  the  following  page),  the  electrical  center  of 
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the  coil,  is  grounded,  a  ground  plana  ZA  ia  sxpactad  in  tba  tabs  as 
a  result  of  the  symmetry  of  tba  coil  and.  the  posh-poll  feed. 

ZIAUsZE M  ItfiflB*  If  the  exciting  frequency  ia  not  too  high. 


i.e. ,  if  its  associated 

wavelength  is  long  com- 

pared  to  the  circuaference 

of  the  coil,  then  the  cur¬ 
rent  at  any  point  of  the 

!  7 

--A 

LX 

coil  at  any  instant  of 

tine  is  constant  and  thus 

the  current  i(t)  is 

independent  of  position. 

The  nagnetic  field  at  1 

position  r  may  be  found 

rig.  i 

/ 

by  use  of  Jnpere's  Law, 

Top  View  of  Exciting  Coil 

(D 


But  by  using  the  law  of  cosines  In  Pig.  1 


x'  =  r  +n  -  2ttL cos 8 

(2) 

and 

2  2  2 

(3) 

2ftr 
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Substituting  |q  (2)  into  Iq  (3) 


slnjp  =  ±  £  -  2afteoafl  -  r2 

2Rr' 


(4) 


Substituting  Eqs  (2)  and  (4)  into  Kq  (1)  and  simplifying 


B 


(r) 


(R2  -  i&cose)d9 
_R2  +  r2  -  2iReoo0>^2 


(5) 


Equation  (5)  is  not  intagrabls  in  olosad  fora. 

By  assuming  that  the  wavelength  associated  with  the  exciting  fre¬ 
quency  was  long  with  respect  to  the  coil  dimension,  Sq  (5)  was  derived 
in  which  B  is  a  function  of  r  and  t  only.  Thus  the  magnetic  flux 
lines  are  sjMetric  about  the  axis  and  consequently  the  lines  of  force 
representing  the  electric  field  E  will  font  closed  concentric  loops 
about  the  axis. 

Faraday's  Law  states  that 

|E  •  dl  =  -  Jafi  •  dl  (6) 

"  8t 

Although  Eq  (5)  cannot  be  integrated  in  closed  font,  nsy  be 

evaluated  at  r  =  0,  which  is  a  point  of  considerable  interest,  yielding 

»«» -  ^  <7> 
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Sine*  B(r,t)  i»  flnlt*  at  r  =  0  and  la  a  continuous  function 
wrtrjrUwN  lnald*  the  coil,  it  will  b*  asaraaed  that  In  the  Interval 
froa  r  =  0  to  r  =  Ar  that  ia  constant  at  any  tlas  t  and 

la  given  fay  Kq  (7) .  Solving  Eq  (6)  for  I  yields 

O  •  r®  *  f-  (8) 

*  •  2«  =  ^  XT2  (9) 


tvzikiim  at 

as  at 


(10) 


Therefore, 


lin  B  =  0 
Ar  ■*  0 


(ID 


Froa  the  preceding  discussion,  it  ia  seen  that  the  probe  aust  be 
positioned  along  the  axis  of  the  tube  to  avoid  RF  pick-up  by  the  probe 
and  to  prevent  the  electrons  froa  being  accelerated  to  the  probe  fay 
the  electric  field. 
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III.  General  Probe  Theory 

Langmuir  wee  the  first  to  do  an  adequate  study  of  the  theory  of 
electron  collection  by  probes.  Later  he  and  Mott-Smith  published  a 
paper  in  October,  1926  (Ref  13)  vhich  is  a  valuable  reference  for  probe 
theory.  Druyvesteyn,  however,  was  the  first  to  show  that  the  second 
derivative  of  probe  current  with  respect  to  probe  voltage  in  the  elec¬ 
tron  retarding  region  is  proportional  to  the  energy  distribution  of 
the  electrons  in  the  plasma.  In  1961  Hedicus  developed  the  theory  of 
spherical  probes  (Ref  11)  using  an  impact  parameter  approach  to  the 
problem.  His  theory  has  the  advantage  over  Druyvesteyn' s  in  that 
although  it  yields  the  same  results,  it  is  much  simpler  to  comprehend. 
To  apply  Medicus'  theory  in  its  entirety  the  data,  however,  obtained 
from  probing  a  plasma  must  come  from  a  small  spherical  probe. 

The  electrons  in  the  plasma  may  assume  any  type  of  distribution 
in  Medicus'  theory,  and  his  theory  will  apply  to  RF  excited  plasmas  as 
well  as  DC  excited  plasmas,  providing  there  is  no  RF  potential  between 
the  probe  and  the  plasma  and  the  probe  is  in  a  field-free  region.  In 
both  the  DC  and  RF  cases,  the  probe  must  be  small  (diameter  on  the 
order  of  0.5  mm)  to  insure  that  the  probe  does  not  distort  the  plasma 
appreciably.  Since  Medicus'  analysis  of  electron  energy  distributions 
depends  on  data  from  probes,  an  understanding  of  the  current  collecting 
mechanism  of  probes  is  Important. 
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Xbt  Probe  Negative  with  Haapect  to  Q|a§  Potential 


Whan  a  probe  is  placad  in  a  plasma  at  a  potantial  strongly  nega- 
tive  with  respaot  to  tba  plasma,  tha  proba  will  rapal  alactrons  and 
attract  ions.  A  positive  ion  sheath  will  build  up  around  the  probe. 
Beyond  the  positive  ion  sheath  there  exists  a  pre-sheath  which  con¬ 
tains  both  positive  and  negative  carriers  and  beyond  this  pre-sheath 
lies  the  undisturbed  plasma.  The  probe  current  at  the  strongly  nega¬ 
tive  potential  will  consist  entirely  of  an  ion  current  as  region  A, 

Fig.  2  indicates. 

As  the  potential 
of  the  probe  is 
made  more  positive 
with  respect  to  the 
plasma  (although  it 
still  remains  nega¬ 
tive  with  respect 
to  the  plasma)  some 
of  the  higher  energy 
electrons  will  over¬ 
came  the  probe- to- 
plasma  retarding 
potential  and  will 
strike  the  probe 

and  be  collected.  In  addition  positive  ions  will  still  be  accelerated 
to  the  probe  so  the  resulting  current  will  be  the  ion  current  minus  the 
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electron  current.  Thus  the  current  collected  by  the  probe  fells  off  es 
region  B,  Fig.  2  indicates.  The  sheath  around  the  probe  will  no  longer 
be  a  pure  positive  ion  sheath  since  electrons  can  penetrate  it,  but 
will  have  the  characteristics  of  a  pre-sheath.  As  the  probe  is  Bade 
even  more  positive  with  respect  to  the  plaaaa  a  point  will  be  reached, 
point  1,  Fig.  2,  in  which  the  rate  of  collection  of  ions  and  electrons 
is  the  sane.  Therefore,  the  probe  will  draw  aero  current.  This  is 
the  potential  at  which  an  Insulated  probe  will  stabilise  when  placed 
in  a  plasma  and  therefore  is  called  floating  potential.  As  was  men¬ 
tioned  in  Chapter  II,  an  ion  will  never  achieve  the  energy  of  an 
electron  because  of  the  comparability  of  ion  and  neutral  gas  atom 
mass;  hence,  the  insulated  probe  will  stabilize  at  a  negative  poten¬ 
tial  because  the  electron  velocity  on  the  average  must  be  greater 
than  the  positive  ion  velocity.  Thus,  more  electrons  will  diffuse  to 
the  probe  per  unit  time  than  positive  ions.  Since  the  probe  is 
insulated,  no  current  can  flow.  The  potential  of  the  probe  will  then 
build  up  to  that  negative  value  at  which  sufficient  mabers  of  the 
electrons  are  repelled  to  cause  the  electron  and  ion  collection  rates 
to  be  the  same.  Theoretically,  the  use  of  the  insulated  probe  was 
one  of  the  earliest  methods  for  determining  plasma  potential  as  it 
was  erroneously  assumed  that  an  insulated  probe  would  assume  the 
potential  of  the  plasma  (Ref  12t315). 

When  the  probe  is  made  more  positive  than  floating  potential, 
more  of  the  electrons  in  the  high  range  of  the  electron  energy  spec¬ 
trum  will  overcome  the  probe- to-plasma  retarding  potential  and  will 
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be  collected.  The  probe  current  will  then  be  equal  to  the  electron 
current  minus  the  positive  ion  current  and  thus  will  appear  as  an 
electron  current  as  region  C,  Fig.  2,  indicates. 

Ihft  Probe  at  flQd  BlgbflC 

Finally,  if  the  probe  is  made  even  more  positive,  it  will  reach 
a  potential  equal  to  plasma  potential,  point  2,  Fig.  2,  and  theoreti¬ 
cally  electrons  and  ions  of  all  energies  will  reach  the  probe,  for 
the  probe  curve  of  Fig.  2,  plasma  potential  is  reached  at  the  point 
where  the  probe  curve  Inflects  as  will  be  shown  in  Chapter  IV.  In  the 
region  ranging  from  the  probe  at  the  potential  at  which  the  highest 
energy  electron  may  be  collected  to  the  probe  at  slightly  less  than 
plasma  potential,  i.e. ,  regions B  and  C,  the  sheath  around  the  probe 
contains  both  positive  and  negative  carriers.  This  region  is  generally 
known  as  the  electron-retarding  region.  At  plasma  potential  there  is 
no  sheath  around  the  probe  because  the  probe  neither  attracts  nor 
repels  electrons  or  ions.  Thus,  the  electrons  and  ions  that  reach 
the  probe  reach  the  probe  by  diffusion  only. 

Making  the  probe  more  positive  than  plfrm  potential  will  cause 
electrons  to  be  accelerated  and  the  positive  ions  to  be  repelled  as 
region  D,  Fig.  2,  indicates.  The  sheath  around  the  probe  will  be  an 
electron  sheath  containing  predominantly  electrons  but  also  those  ions 
with  energies  high  enough  to  overcame  the  probe- to-plmsma  potential. 
Increasing  the  potential  of  the  probe  sufficiently,  region  K,  Fig.  2, 
causes  all  ions  to  be  repelled  and  also  oauses  an  exclusively  electron 
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sheath  to  font  around  ttaa  proba.  A  surra  such  aa  that  in  Fig.  2  nay 
ba  automatically  plottad  aa  will  be  diaeuaaad  in  Chapter  T. 

Region  of  Interest  for  Iggggy  nistrifamtion 

The  primary  region  of  interest  for  determining  energy  distribu¬ 
tions  is  the  electron  retarding  region.  As  has  been  shown,  the  magni- 
tude  of  the  probe  current  varies  as  the  probe- to-pla saa  potential  is 
varied.  This  change  in  the  electron  retarding  region,  as  already  dis¬ 
cussed,  is  a  result  of  the  fact  that  electrons  with  lower  energies  are 
able  to  overcome  the  probe- to-pla  ana  potential  as  the  potential  is  made 
■ore  positive,  until  finally  at  plasma  potential  electrons  of  all 
energies  will  be  collected.  Intuitively  then  there  should  be  a 
relationship  between  the  electron  energy  distribution  of  the  plasma, 
the  probe  current  and  the  probe  potential.  The  next  chapter  will  show 
how  they  are  related. 
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17.  sobsslSAl  ttatu  zbtficx 

Medicus'  spherical  probe  theory,  which  he  proposed  in  the  Journal 
of  implied  Phralca  in  December,  1961 ,  was  used  exclusively  in  this 
project  and  will  be  discussed  in  this  chapter.  This  chapter  parallels 
his  spherical  probe  theory  (Ref  11)  in  shoving  the  relationship  between 
the  energy  distribution  of  the  electrons  in  the  plasma,  the  probe  cur¬ 
rent  and  the  probe  potential,  but  this  author  has  found  it  more  advan¬ 
tageous  to  this  paper  to  derive  all  equations  in  terms  of  energy 
distribution  functions  rather  than  velocity  distribution  functions 
as  was  done  by  Medicus.  The  equation  used  by  Medicus  for  determining 
current  density  will  be  derived  by  this  author  on  the  basis  of  electron 
scattering.  The  location  of  plasma  potential  on  the  probe  curve  will 
be  shown  as  a  logical  extension  of  Medicus*  spherical  probe  theory. 

Mott-Smith  and  Langmuir  in  1926  devised  equations  for  electron 
collection  in  an  article  In  Physical  Review  but  suggested  an  alternate 
attack  which  consisted  of  subdividing  the  flow  of  electrons  into 
parallel  moving  swarms  of  electrons  (Ref  1 3i751 ) .  Medicus  has  used 
this  approach  in  developing  his  theory  of  electron  collection  by  the 
spherical  probe. 

Probe  Collection  flf  aq  Isotropic  Distribution  flf  MPBflflaflTgflUC  ElfStfPPa 

If  a  probe  having  a  radius  r^  is  immersed  in  a  hypothetical 
plasma  containing  mono  energetic  electrons  and  is  at  the  potential  of 
the  plasma  in  an  actual  case  it  will  collect  ions  and  electrons  which 
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diffuse  to  the  probe.  The  foot  that  ions  diffuse  to  tha  probe  whan  it 
is  at  plasne  potential  and  are  attraoted  to  the  probe  when  it  is  nega¬ 
tive  with  respeet  to  the  pi  sms  will  be  ignored  in  the  following 
derivations  as  the  effect  of  the  positive  ion  current  will  later  be 
quantitatively  shown  to  be  null,  for  the  purpose  of  this  derivation 
then  the  positive  ions  will  be  aasuaed  to  be  in  the  plasma  for  the 
sole  purpose  of  poreserving  space  charge  neutrality  and  will  not  be 
collected  by  the  probe. 

The  current  collected  by  the  probe  will  therefore  be  purely  an 
electron  current  and  will  be  equal  to  the  product  of  the  collecting 
area  of  the  probe  and  the  electron  current  density.  If  it  is  aasuaed 
that  the  distributions  of  the  electrons  in  the  hypothetical  planes  are 
isotropic  and  eaoh  electron  has  the  aaaa  speed  v  then  the  electron 
current  density,  as  will  be  shown  in  the  paragraph  that  follows,  ia 
9BX  where  q  is  the  electronic 
charge  end  n  is  the  electron 
number  density.  Tbs  collecting 

area  of  the  probe  is  the  surface 

2 

area  of  the  probe,  i.e.,  4xr^ 

as  fig.  3  suggests. 

By  an  electron  scattering 

approach,  tha  currant  density 

ia  quite  easily  determined.  In  fig.  4  (shown  on  the  following  page)  S 

is  a  unit  of  area  normal  to  the  Z  direction.  If  2  represents  tha 

8 

macroscopic  electron  elastic  scattering  cross  section  of  the  neutral  gas 


fig.  3 

Currant  Collection  in  an 
- Isotropic  Plasne 
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•teas*  than  the  maker  of  alartmi  elaatlaaUy  mtUwd  oat  of  volane 
elonent  r2sined9dfdr  la  anrsgr2alaed»dfdr.  Mot  all  of  theao  eleo- 
trona  will  roach 
aroa  8  aloeo 
km  will  bo 
attooaatod  on 
route  bpr  olaatlo 
aeattorlng  and 
alao  bacanaa  tha 
aeattorlng  ia 
aaeunod  to  bo 
iaotropic.  Tha 
mihor  aoattorod 
out  of  tho  volme 
olanant  and 
reaching  S  la 

n^rV^a* 

thua  -  .r.  —  alaOeoafldModr.  Tha  fact  that  eaae  olootrona  eoold  bo 
4ar2 

iaalaatioally  aoattorod  ia  lgnorod  bacanaa  for  no  at  Talnos  that  ▼  will 
aeeune,  tho  aaorosooplc  lnalaatie  aeattorlng  eroaa  aoction  ia  anall. 

Tha  currant  doneitj  J  atr iking  S  fron  tho  poaitive  Z  direction  ia 
thoroforo  given  by 


4 


Pig.  4 

Area  S  into  tdiich  Eloc trona 
arc  Soattorod 


1 '  £  Jf  J?  “  Vv 


(12) 
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that  tha  area  3  is  far  from  the  boundaries  of  the 
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If  it  is  assaned 

plana,  the  upper  Unit  of  r  nay  be  set  at  infinity.  Therefore 

J  =  C  £ Jo  ^  V’*8*  co sealnedOd^dr  (13) 


and  upon  integration 


J  =  3EI 
4 


(14) 


Since  I,  the  probe  current  is  the  product  of  the  collecting  area  and 
the  current  density 


j  _ 

4 


(15) 


or 


I  »  XX  qpr 
P 


(16) 


Ezaht  flaUjatto  aL  a  EMaliil  Smub  al  vaacamtMU£  neetrona 

If  the  s an  spherical  probe  is  placed  in  a  nono-directional 
parallel  swam  of  nonoeoergetio  electrons,  Fig.  5,  and  the  probe  poten¬ 
tial  is  the  sans  as  plasna  poten¬ 
tial  the  current  collected  by  the 
probe  is  again  equal  to  the  prod¬ 
uct  of  the  collecting  area  and 
the  current  density  which  are 
now  different  fron  the 


3  0 

Fig.  5 

Current  Collection  in  a 
Parallel  Swam  of  glectrons 
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isotropic  cass.  Tbs  collecting  ares  for  tbs  paroUsl  swarm  is  tbs 

cross- sectional  area  of  tbs  probe  (fig.  5  on  tbs  prsesding  page) , 

2 

i.e. ,  xTp.  Tbs  current  density  is  qnv.  Thus, 

1  =  xr^inv  (17) 

This  is  tbs  sane  equation  as  Sq  (16).  Henoe,  a  spherical  probe  in  a 
plasm  of  monoenergetio  electrons  can  be  analysed  as  a  sphsrioal  probe 
in  a  Bonoenergetic  parallel  swarm  of  eleotrons.  By  further  extension 
of  this  idea  a  sphsrioal  probe  in  any  plasm  can  be  analysed  as  a 
spherical  probe  in  many  parallel  swarm  of  eleotrons  in  which  each 
parallel  swarm  corresponds  to  a  particular  energy  range. 

iff sot  Probe  Potential  on  Bffsctlve  Collecting  Area 

If  a  potential  is  applied  to  the  probe  that  is  negative  with 
respect  to  plasm  potential,  the  electrons  will  be  repelled  from  the 
probe.  If  a  positive  potential  is  applied  with  respect  to  the  plasm 
potential,  the  electrons  will  be  attracted  to  the  probe.  Their 
behavior,  when  in  the  presence  of  the  probe,  can  best  be  examined  by 
considering  a  monoenergetic  parallel  swarm  of  electrons  and  by  deter¬ 
mining  their  Impact  parameters.  The  Impact  parameter  of  an  electron 
in  the  presence  of  the  probe  is  the  distance  by  which  the  electron 
would  miss  the  center  of  the  probe  if  there  were  no  Coulomb  force 
between  them.  The  impact  parameter  for  graslng  incidence  on  the 
probe  is  the  Impact  parameter  of  an  electron  which  just  touches  the 
surface  of  the  probe  as  it  passes  by  the  probe. 
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If  Pg  reprs— nts  the  lnpaet  parameter  for  graslng  incidence  on 
the  probe,  then  electrons  with  lnpaet  parameter  lees  than  p^  will  be 
collected  by  the  probe  end  eleotrons  with  Inpact  perenetere  greeter 
than  will  not  be  oolleoted.  Figure  6  illustrates  two  situations. 
Fig.  6a  is  the  case  of  electron  acceleration,  i.e. ,  V  >  0  where  V 
is  the  poten¬ 
tial  of  the 
probe  refer¬ 
enced  to 
plane  poten¬ 
tial.  Figure 
6b  is  the 
oase  of  elec¬ 
tron  retarda¬ 
tion,  i.e., 

V  <  0.  the 
different 
paths  the 

eleotrons  assume  in  the  two  cases  is  due,  of  course,  to  the  Coulomb 
force  between  the  electrons  and  the  probe.  In  a  parallel  swam  of 
electrons,  if  the  probe  is  at  a  potential  other  than  plasm  potential, 
the  collecting  area  will  not  be  the  arose- sectional  area  of  the 

2 

spherical  probe  but  will  be  a  function  of  p  specifically  *p  .  As 

D  6 

Fig.  6  Indicates,  in  the  case  of  a  probe  potential  negative  with  respect 
to  plasm  potential,  the  collecting  area  is  less  than  the  cross-sectional 


v>o 

v<o 

(•) 

(b) 

Inpact  Para) 

Fig.  6 

leters  for  erasing  Incidence 

(Fron  Ref  11*2312) 
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irn  of  the  probe  since  p^  <  rp.  In  the  out  of  tho  probo  at  a  poten¬ 
tial  positive  with  respect  to  the  pi  sms,  the  collecting  area  is  greater 
since  pg  >  rp.  in  equation  for  probe  current  in  a  nonoenergetlo  paral¬ 
lel  swan  of  electrons  for  a  spherical  probe  at  any  potential  nay  again 
be  written  as  the  product  of  the  collecting  area  and  current  density. 

The  collecting  area  of  the  probe  for  probe  potentials  other  than  plsmsa 

2 

potentials  is  xp^  within  limitations  which  will  be  discussed  on  page 
25.  The  current  density  is  again  qnv.  Therefore, 

I  =  *Pgqnv  (18) 

which,  following  the  reasoning  used  previously,  is  not  only  the  equation 
for  current  collection  in  a  nonoenergetlo  parallel  swan  of  electrons, 
but  is  also  the  equation  for  current  collection  in  a  plana  of  nono- 
energetic  electrons. 

p  nay  be  calculated  as  a  function  of  probe  radius,  probe  poten- 

o 

tial,  and  electron  energy  by  considering  the  equations  of  conservation 
of  energy  and  conservation  of  angular  nouentun  about  the  probe  and  by 
assuming  the  potential  drop  between  probe  and  plaaaa  occurs  across  the 
sheath. 

By  the  law  of  conservation  of  energy,  the  kinetic  energy  of  an 
electron  at  the  probe  equals  the  sun  of  the  kinetic  energy  of  the 
electron  at  the  sheath  and  its  potential  energy,  l.e., 

2  *vp  a  \  wr\  +  qV  (19) 
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where  vp  is  the  velocity  of  the  electron  at  probe; 

▼  is  the  velocity  of  tha  alactron  at  ahaath. 

▼p,  however,  is  composed  of  radial  and  tangential  oonponanta  to 
tha  probe,  i.a. , 


2  2  2 

▼„=▼'+  vf, 

p  rp  tp 


(20) 


vhara  is  the  radial  ooaponant  of  velocity  at  proba; 

Ttp  is  tha  tangantial  ooaponant  of  velocity  at  proba. 
v^,  however,  aust  equal  aero  in  order  for  an  alactron  to  grass 
tha  proba.  Thus, 


(21) 


Substituting  Eq  (21)  into  Xq  (19) 

£  ■'tp  =  2  +  «V  (22) 

Since  angular  aoaantua  about  tha  proba  mist  be  conserved,  tha 
angular  noaantm  of  tha  alactron  at  tha  sheath  aust  equal  tha  angular 
■onantun  at  tha  proba. 


r  x  nv  *  r  x  mr  (23) 

s  p  p 

For  grasing  ineidanoa  at  tha  proba, 

Pg*%  =  V^tp  <24> 
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and  thus, 


P_T 

\p  =  ^ 

'P 


Substituting  lq  (25)  into  Sq  (22)  yields 


1,  £gli  _  1_2 


2  _2  2 
'p 


Dividing  lq  (26)  by  1  — £  yields 

rp 


2  2  2rfqV 

p  =  r_  +  r 

*  p  «rj 


or 


(25) 


(26) 


(27) 


(28) 


where  1  is  the  energy  of  the  electron  et  the  sheath  edge  which  is  the 
energy  of  the  electron  in  the  undisturbed  plana.  Equation  (28)  nay 
also  be  written  as 


r>  *  §> 


(29) 


where  U  is  the  kinetic  energy  of  the  electron  in  electron  volts. 
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£g  IK  Potltifi  jbA  MfatAy  Prtht 

It  this  time  it  is  worthwhile  to  st— ins  Sq  (28)  in  toms  of 
Fig.  6.  When  tbs  potential  of  tbs  probe  is  sore  positive  than  pi  a— 
potential,  pg  is  greater  than  r^  for  all  energies.  Bqoation  (28) 
demonstrates  that  Pg  at  a  fixed  potential  approaches  inf  laity  for 
energies  approaching  aero.  This  is  not  logical  except  in  the  case  of 
the  probe  with  an  Infinite  sheath.  Since  the  potential  drop  between 
probe  and  planes  is  assumed  to  occur  across  the  sheath,  only  those 
eleotrona  reaching  the  sheath,  i.e. ,  those  with  Pg  *  rg,  are  attracted 
to  the  probe.  Therefore  in  the  elsctrom-aocsleratlng  region  the  upper 
Unit  of  Pg  equals  rg.  Thus  for  0  <  S  <  1^,  pg  must  equal  rg. 
Substituting  Pg  a  rg  into  Sq  (28)  yields 


*J“rJ(i  +  ai) 


(30) 


Therefore, 


The  region  of  interest  in  deternlnlng  electron  energy  distributions, 
however,  is  for  probe  potentials  less  than  plana  potential. 

Uhan  the  potential  of  the  probe  is  less  than  plasma  potential  Pg 
is  indeed  less  than  r^  far  |qV|  £  K  <  •  as  Sq  (28)  indicates.  Thus 
the  collecting  area  is  less  than  the  cross- sectional  area.  The  Units 
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are  sat  on  X  sine*  an  infinitely  faat  alaetron  has  a  gracing  Inpact 
par— star  p  =  r  while  —  alaetron  with  energy  E  3  |qV|  has  a 

O  m 

gracing  lapaet  par— ter  p(  =  0.  Klaotrona  with  1  <  |qV|  will  not 
be  eolleetad  by  the  probe. 

QgCigg&Jfig  QtpgggX  Cnrrcnt  Equation  jg  Eleatnan  gg^gQ^gg  Region 

By  aabatitutlng  Eq  (28)  into  Bq  (18),  Sq  (18)  my  now  be  written  as 

1  =  *r2(1  +  Stf)qnv  (32) 

where  V  will  be  restricted  to  ▼slues  negative  with  respect  to  plasm 
potential.  Converting  Eq  (32)  to  a  function  of  energy, 

I  =  crpqn(J)1/2(1  +  ^Jl1/2  for  V  <  0  (33) 

General  enargy  distributions  nay  now  be  considered  by  superposition 
of  lnfinitesiaal  nonoenergetie  contributions.  If  F(I)  is  defined  as 
the  fraction  of  electrons  par  unit  energy  Interval,  than  F(S)dS  is  the 
fraction  of  electrons  with  energies  between  E  and  X  +  dE  and 
nF(E)dE  is  the  uuaber  density  with  energies  between  E  and  E  +  dE. 

The  contribution  to  the  probe  currant  for  an  energy  range  free  E  to 
E  +  dE  is  therefore 

dl  =  *r2qn(J)1/2(1  ♦  f^X^dJdE  for  V  <  0  (34) 

Since  Eq  (34)  is  restricted  to  probe  potentials  negative  with  respect 
to  plasm  potential  the  total  current  at  a  particular  potential  nay  be 
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found  tgr  integrating  lq  (34)  over  *11  energies  of  •loetrons  that  eon  bo 

oollootod  hf  tbo  probe.  For  •  particular  V  <  0  tho  lowest  energy 

electron  that  iqr  bo  oollootod  ia  By>  where  By  approaches  qV  as  a 

liait,  i.e.,  11b  (By  4-  AS)  *  qV.  Thus 

AS  **  0 


I  * 


C(1  +  ^)I1/2F(I)dB 

*7  S 


where  K  =  nr' 


(35) 


of  Electron  laargy  niirtrtfcuttaa  ranotton 
By  differentiating  lq  (35)  twice,  it  will  now  bo  shown  that  tho 
distribution  function  can  bo  expressed  in  tans  of  tho  second  deriva¬ 
tive  of  probo  current  with  respect  to  probe  voltage. 


S  =  *D|)*1/2F<*>«  06) 

nr 

dV  =  (37) 


Sines  By  approaches  qV  as  a  Unit 


dly  =  qdV 


(38) 
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AL  *  1  41 
dly  q  dv 


Substituting  Sq 


dS 


q  d«y 


(38)  into  Sq  (40)  yields 


d^  « 


(») 

(40) 


(41) 


and 


^1  =  1  & 
dsj  q2  dv2 


(42) 


Substituting  Sq  (39)  into  Sq  (37)  yields 


(43) 


Differentiating  Sq  (43)  with  respect  to  S^  and  simplifying  yields 


=  K(Sy)“l/2F(By) 

dE„ 


(44) 
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Substituting  Eq  (42)  into  Eq  (44)  yields 


&  =  Kq2(V"l/2F(V  (49) 

ar 

or  P(I  )  =  K’  &  (*y)1/2  (46) 

V  dV2 

where  K'  =  Kq2. 

Tbs  significance  of  Sq  (46)  is  that  tbs  fraction  of  slsetrons 
having  a  particular  energy  Ey  (whera  By  approaches  qV  as  a  llolt) 
per  unit  energy  Interval  can  be  detenained  by  Multiplying  the  second 
derivative  of  probe  current  with  respect  to  probe  potential  by  the 
square  root  of  the  energy  Ey.  Ey  and  the  second  derivative  of  probe 
current  can  be  obtained  from  the  second  derivative  of  a  probe  curve 
automatically  plotted  from  data  taken  by  a  spherical  probe.  If  Eq  (46) 
is  applied  throughout  the  electron-retarding  region  th#  energy  dis¬ 
tribution  of  the  electrons  may  be  detenained  as  will  be  discussed  More 
fully  in  Chapter  TI. 

location  £laaaa  Potential 

The  fact  that  plasma  potential  occurs  at  the  inflection  point  of 
a  probe  curve  will  now  be  shown,  is  has  been  Mentioned  on  page  26, 
when  the  probe  is  at  a  positive  potential,  those  electrons  whose  ener¬ 
gies  are  less  than  E^  will  have  an  impact  parameter  for  grazing 

incidence  of  p  =  r  .  For  an  infinitely  thin  sheath,  i.e.,  when 
a  8 
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r#  ■*  rp,  1^  •  as  Iq  (31)  indicates.  Banes  An  r#  -*  rp  alactrona 

of  all  energies  will  have  p_  =  r  and  the  proba  current  for  7  >  0 

a  ■ 

la  given  by 


I  *  J*«r^qn(|)1/2l1/2F(»)d*  for  V  >  0  (47) 

pa  '*  rp 

2 

ainca  tba  collecting  area  is  now  xr#. 

Tha  above  aquation  is  indapandant  of  V)  banoa  tha  currant  will 
aaturata  above  plana  potential  for  tha  oaaa  of  tha  infinitely  thin 
sheath.  Since  r  -»  r  ,  tha  value  at  which  tha  current  will  saturate 

8  p 

is  tha  value  of  tha  current  at  plana  potential  because  at  plaana 
potential  1^  =  0  and  V  =  0  in  lq  (24) . 

For  the  other  extraae  ease  of  an  Infinitely  large  sheath  •*  0 
in  lq  (31).  8inee  rfl  is  infinitely  large,  lq  (28)  gives  the  correct 
value  of  Pg  for  all  7.  Thus  the  currant  is  now  given  by 


I  =  J**rpqn(J)1/,2(1  +  ^) (l1/2)P(l)dl  for  V  >  0  (48) 


Upon  integration,  I  will  be  in  tha  fora  of  I  =  a  +  bV  where  a  and 
b  are  oonstanta  (Ref  11s2515).  Since  lq  (48)  at  7=0  is  the  sons 
as  lq  (35)  at  7=0  and  lq  (48)  is  linear  in  7,  the  portion  of  the 
proba  curve  in  the  accelerating  region  is  tangent  to  tha  probe  curve 
in  the  retarding  region  and  is  linear. 
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Actual  probes  hava  neither  infinitely  thin  sheaths  or  infinitely 
large  sheaths.  Therefore  actual  probe  curves  in  the  electron- 
accelerating  region  Bust  lie  between  the  two  limiting  cases,  Fig.  7. 
Since  the  upper  limit  of  the  probe  curve  for  V  >  0  is  linear  and 
tangent  to  the  curve  for  V  <  0  the  actual  probe  curve  at  most  must 
be  linear.  It  more  probably  will  be  concave  downward,  however. 

The  probe  curve  in  the  electron-retarding  region  will  now  be 
shown  to  be  concave  upward.  For  the  monoenergetic  case,  it  has  been 
shown  that 


I  =  nvq(1  +  |r)*rp 


(32) 


Thus  for  a  constant  E,  I 
varies  linearly  with  V.  At 
the  potential  V  where 
E  =  -qV  the  current  will 
be  zero.  When  V  =  0,  1 

is  at  its  maximum  in  the 
electron  retarding  region, 
i.e., 

2 

I  =  nvqurp  (16) 


V  — » 


rig.  7 

Actual  Probe  Curve  with 
Two  Halting  Cases 


If  one  considers  a 

hypothetical  plasma  containing  electrons  of  three  energy  groups  , 
E2,  and  E3  in  which  (a)  E1  >  E2  >  E3,  (b)  the  number  density  of 
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each  group  is  ,  n^,  and  n^,  respectively,  and  (c)  the  velocity 

of  eaoh  group,  i.e.,  V1 ,  V2»  and  V  respectively,  then  a  probe 
curve  such  as  the  one  in  Fig.  8  may  be  constructed  with  the  use  of 
Bq  (32)  •  This  idea  may  be  extended  to  a  plasma  containing  all  energies 
with  the  resulting 
figure  concave 
upward.  Thus 
plasma  potential 
does  indeed  cor¬ 
respond  to  the 
inflection  point 
in  a  probe  curve 
since  it  marks  the 
transition  from  a 
curve  that  is 
concave  upward  to 

one  that  is  either  linear  or  more  probably  concave  downward.  Plasma 
potential  also  must  correspond  to  the  point  where  the  second  derivative 
of  the  probe  curve  equals  zero  since  the  second  derivative  of  any  curve 
at  its  inflection  point  is  zero. 

Calculation  si  Plasm  Penalty 

From  Fig*  6  it  can  be  seen  that  the  current  at  plasma  potential 
for  a  plasma  having  electrons  of  energies  E^,  E2>  and  is 
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given  by 

3 

I  =  *rpq  I  Vi  (49) 

i=1 

Therefor*!  for  a  plasma  having  N  energy  groups 

N 

1  =  "pq  I  Vi  <50> 

i=1 


But  since 

N 

I  Vi 

V  S  tel  .  ■■  (51 ) 

N 

E"i 

i=1 

N 

and  n  =  L  ni  (52) 

i=1 

N 

and  M  =  £  n^Vj  (53) 

i=1 


I  =  itr2qnV  (54) 

2  — 

or  n  »  i/xr^qv  (55) 

Therefor* ,  if  the  radius  of  the  probe  and  the  current  at  plasma  poten¬ 
tial  are  known,  plasma  density  may  be  calculated  by  first  coeqputing  the 
average  electron  velocity  and  then  using  Eq  (55) . 
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V. 

The  following  stops  wore  undortokon  tjr  tbs  author  In  preparing  tho 
equipment  for  data  taking: 

(1)  Design,  construction  and  mounting  of  tho  plasma 
tube; 

(2)  "Bake-out"  of  tho  tuba  and  subsequent  filling 
with  argon; 

(3)  Installation  of  the  plasm*  exciting  equipment) 

(4)  Installation  of  tho  measuring  equipment; 

(5)  Positioning  of  tho  probe  along  tho  axis  of  tho 
•setting  ooil; 

(6)  Cleaning  of  tho  probe. 

£lgggg  Tube 

Tbe  tube  that  contained  tho  plasma  is  shown  in  Pig.  9  on  the  fol¬ 
lowing  page.  The  tube  was  11.5  cm  in  diameter  and  20.5  cm  high.  It 
was  designed  with  a  champagne  bottle  shape  on  both  ends  to  prevent 
collapsing  at  high  vacuums.  The  tube  and  the  entire  system  were  tested 
with  a  helium  leak  detector  at  a  vacuum  of  10"^  torr  and  no  leaks  were 
found. 

Since  there  is  no  physical  way  in  which  to  reference  a  probe  to  a 
plasma  directly,  a  reference  electrode  must  be  located  in  the  plasma 
for  all  probe  measurements,  for  a  DC  plasma  the  probe  may  be  refer¬ 
enced  to  one  of  the  electrodes  already  in  the  tube,  i.e. ,  the  anode  or 


34 


GNE/Phys/63-3 


Rtftrnu 


Pig.  9 

Plan*  Tub* 
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oathode.  For  an  electrodeless  discharge  (EF  excited  plasma)  a  reference 
electrode  sust  be  placed  in  the  tube. 

The  reference  electrode,  Fig.  9,  used  in  this  experiment  was  a 
1 0-cm- In-dlameter  circular  Monel  disk  0.5  an  thick.  Monel  was  chosen 
because  it  is  easy  to  work  with  and  is  non-magnetic.  The  rather  large 
size  of  the  reference  electrode  was  selected  for  two  reasons.  First, 
it  was  easier  to  construct  than  a  smaller  one  would  be.  Secondly,  Fetz 
and  ueohsner  (Ref  4t251 )  reccanended  in  their  paper  on  the  determination 
of  plasma  potential  of  an  RF  excited  plasma  that  the  ratio  of  surface 
area  of  the  reference  electrode  to  the  surface  area  of  the  probe  be  at 
least  lOOOtl.  The  ratio  in  this  case  was  22,0000. 

The  physical  reason  for  the  large  ratio  between  probe  and  reference 
electrode  area  is  obvious.  When  the  probe  draws  electrons  from  a  piaster, 
the  reference  electrode  draws  an  equal  number  of  positive  ions.  If  the 
reference  electrode  were  the  same  size  as  the  probe,  when  the  potential 
between  probe  and  plasma  is  made  more  positive,  the  probe  would  only 
draw  an  increasing  amount  of  electron  current  until  ion  saturation  of 
the  reference  electrode  occurred.  When  ion  saturation  of  the  reference 
electrode  is  reached,  electron  saturation  of  the  probe  must  be  reached 
because  the  probe  cannot  remove  more  electrons  from  the  plasma  than  the 
reference  electrode  can  remove  ions.  In  examining  the  typical  probe 
curve  in  Fig.  2,  it  is  evident  that  ion  saturation  of  the  reference 
electrode  would  distort  the  shape  of  the  electron  current  of  the  probe 
drastically.  If  the  reference  electrode  is  made  many  times  larger  than 
the  probe,  the  reference  electrode  will  then  increase  its  ion  collecting 
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capability  by  a  factor  approximating  tha  ratio  of  tha  rafarance  alac- 
troda  area  to  tba  araa  of  tha  proba.  If  tbla  it  dona,  tha  rafaranoa 
alactrode  will  not  saturate  bafora  alaotrcn  saturation  of  tha  proba 
occurs  and  hanca  will  hava  no  effect  on  tha  electron  drawing  capability 
of  tha  proba. 

A  1 -cm- in-diameter  hole  was  cut  in  tha  center  of  tha  rafarance 
alactrode  to  allow  tha  proba  to  pass  through.  Tha  two  rows  of  seven 
0.1~em-wlde  evenly  spaced  slits  in  tha  rafarance  electrode  (Tig.  9) 
ware  cut  for  the  purpose  of  preventing  eddy  currants.  Spot  welded  to 
tha  rafarance  electrode  ware  two  0.1-ca  inside  diameter  Monel  sleeves 
7  cm  apart,  into  each  sleeve  was  fitted  a  Kovar  wire,  8  cm  long  and 
0.1  cm  in  diameter.  These  two  wires  suspended  tha  rafarance  electrode 
in  tha  bottle  and  protruded  1  cm  through  tha  glass  tuba  to  provide 
electrical  contacts  so  that  tha  proba  could  be  referenced  to  ground. 

Tha  proba  consisted  of  a  0.3<1  ass- in -disaster  tungsten  sphere 
suspended  on  a  tungsten  wire.  Surrounding  tha  wire  was  a  quartz 
capillary  to  prevent  electrons  and  ions  from  being  collected  by  the 
wire  in  addition  to  being  collected  by  the  probe.  The  probe  was 
designed  to  be  capable  of  being  raised  to  within  2.5  cm  of  the  reference 
electrode  and  lowered  to  9  cm  from  the  reference  electrode. 

The  plassui  tube  was  blown  by  the  glass  blower  of  the  Klectronic 
Technology  Laboratory  and  mounted  on  a  pump  stand.  Figure  10,  shown 
on  the  following  page,  is  a  schematic  showing  the  placement  of  the 
tube  on  the  pump  stand  and  the  valve  arrangement  used  for  filling 
the  tube. 
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X12BUB  ftajm 

The  argon  gas  that  was  used  to  fill  tha  plasma  tuba  us  ad  in  this 
investigation  was  supplied  'ey  tha  1 .1  liter  bottle  shown  below  in 


diameter  and  17  cm  long.  This  thin  Kovar  tuba  was  connected  to  valve  3. 

The  bottle  was  sealed  by  the  manufacturer  by  naans  of  a  glass  tip 
as  seen  at  the  junction  of  the  glass  tubing  and  gas  bottle  in  Fig.  10. 

In  order  to  allow  the  gas  to  escape,  this  tip  had  to  be  broken.  Before 
the  glass  tube  was  sealed  to  the  Kovar  tube,  an  iron  slug  was  placed  in 
the  glass  tube.  This  slug  was  moved  up  against  the  glass  tip  (at  the 
appropriate  tine)  by  means  of  a  magnet  until  the  tip  broke.  To  insure 
a  clean  system,  this  iron  slug  must  be  encapsulated  in  glass  before 
it  is  inserted  into  the  glass  tubing  in  order  to  eliminate  outgassing 
from  the  slug. 
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Valves  1,  2,  and  3  were  designed  by  the  Electronic  Technology 
Laboratory  for  high  vacuum  work.  They  have  the  characteristic  of  being 

a 

able  to  hold  a  vacuum  of  10  torr  and  can  withstand  at  least  120 
"bake-outs". 

To  insure  an  impurity-free  system,  the  tube,  valves,  and  connecting 
pipes  had  to  be  baked  before  the  tube  was  filled  with  argon.  Valves  1, 
2,  and  3  were  opened,  the  pumps  were  put  in  operation,  and  an  electrical 
oven  was  lowered  over  the  equipment.  The  tube  was  baked  for  10  hours 
at  450°C  until  all  outgassing  ceased.  The  oven  was  raised,  all  valves 
were  closed,  and  the  iron  slug  was  then  rammed  up  against  the  glass 
tip,  breaking  it  and  allowing  gas  to  fill  the  glass  tubing.  Valve  3 
was  opened  and  then  closed  again  to  allow  gas  to  fill  the  volume  between 
valve  3  and  2.  Valve  2  was  then  opened,  allowing  gas  to  fill  the  plasma 
tube. 

Inciting  Coll 

The  plasma-exciting  coil  used  in  this  investigation  was  a  one- turn 
copper  coil  16.5  cm  in  diameter.  The  copper  coil  was  made  by  bending  a 
52  x  17.7  x  0.2  cm  sheet  of  copper  in  the  form  of  a  cylinder.  The  ends, 
however,  were  not  completely  closed.  A  gap  of  1 .5  cm  remained  and  the 
two  ends  of  the  copper  coil  were  bent  perpendicular  to  the  coil  surface 
to  create  two  0.75- cm- wide  flaps  as  shown  in  Pig.  11  on  the  following 
page.  This  was  done  so  that  a  variable  mica  capacitor  could  be  placed 
across  the  coil  and  to  provide  a  convenient  place  to  attach  the  power 
input  coils.  One  RG-8/U  shielded  coaxial  cable  was  bolted  and  soldered 
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to  Moh  aid*  of  tho  ooil  at  tbs  slit  to  provide  tbs  power  loads  to  tbs 
coil.  Tbs  tank  circuit 
consisting  of  tbs  varia¬ 
ble  capacitor  and  tbs 
copper  coil  was  tuned 
to  the  transmitter 
frequency,  which  was 
21.3  ac.  This  fre¬ 
quency  was  selected 
beoause  it  lies  in 
the  anatour  band  and 
because  its  associated 
wavelength  is  large 
oonpared  to  the  dimensions  of  tho  coll,  i.e. ,  the  ratio  of  wavelength 
to  coil  circumference  is  about  30  to  1.  Because  of  this  fact  and  the 
fact  that  the  ooil  was  fed  push-pull,  the  current  at  any  point  on  the 
coll  at  any  instant  of  time  was  constant)  hence,  the  field  at  the 
center  of  the  coil  was  aero  and  tbs  ground  plane  mentioned  in  Chapter  H 
existed  in  the  tube. 


Power  Source 

The  transmitter  used  as  a  source  of  power  in  exciting  the  plasma 
waa  a  coneercial  Viking  Ranger  made  by  E.  F.  Johnson  Company.  It  had 
a  single-ended  output  rated  at  65  watts.  The  RF  power  waa  fed  by  an 
RG-8/U  coaxial  cable  to  a  T-splioa  contained  inside  a  grounded  aluminum 
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box.  The  output  cables  frou  the  T-aplioa  vara  tba  ones  bolted  and 
soldered  to  the  coil  in  Fig.  11,  as  Mentioned.  The  lengths  of  the 
coaxial  cables  were  0.65  asters  and  5.30  asters. 

The  physical  length  of  wire  corresponding  to  a  half  wave  length 
of  current  or  voltage  is  given  by 


X  =  -C.  x  VF  (56) 

2v 

where  \  =  physical  wave  length  of  wire 

C  =  velocity  of  light  in  free  space 

V7  =  velocity  factor  which  is  defined  as  the  ratio  of  wavelength 
in  wire  to  wavelength  in  free  space 
v  =  frequency  of  operation. 

For  RG-8/U  the  velocity  factor  is  0.66  (Ref  15t348). 


Hence, 


(57) 


X  =  4.65  mtrs  (58) 

The  two  coaxial  cables  used,  therefore,  differed  in  length  by  one- 
half  wavelength.  Thus  at  any  instant  of  tine  the  current  and  voltage  at 
the  coil  end  of  one  cable  were  180°  out  of  phase  with  the  current  and 
voltage  at  the  coil  end  of  the  other  cable.  By  this  judicious  selection 
of  lengths,  the  single-ended  feed  of  the  transmitter  was  converted  to 
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pub-pull.  This  arrangemm: it  for  converting  fron  single  ended  to  pneb- 
pull  food  ms  adopted  bc~ense  of  expediency.  in other  arrangement 
whieh  will  bo  — Btloaod  in  the  following  paragraph  did  not  work  satis¬ 
factorily. 

in  ettoapt  was  made  to  food  the  elngloi  ended  output  of  the  trans¬ 
mitter  to  the  prlaary  of  a  traneforaer  having  a  center-tapped  secondary 
in  order  to  convert  the  single-ended  food  to  puab-pull.  1  good 
Impedance  natch  could  not  be  node  and  the  roeultlng  pie—  wi  eo  weak 
that  probe  curves  could  not  be  taken,  i  laboratoxy-beilt  generator 
(not  designed  for  this  saqporinont)  with  a  pnsh-pnll  output  was  tried* 
but  it  had  no  IF  ground.  Balanced  eaounte  of  power  ooold  not  be  fed 
froa  this  tranaaltter  into  each  side  of  the  coil  since  the  oeoillator 
was  self-exciting.  In  this  oase  the  plum  was  quite  donee*  but  a 
large  eaount  of  BF  interference  was  picked  up  by  the  probe  so  that 
the  oondition  of  no  BF  potential  between  probe  end  plaaaa  oould  not 
be  net}  hence*  no  aeaalngful  probe  curves  oould  be  taken.  The  anther 
started  building  a  crystal  controlled  oscillator  having  a  pneh-puH 
power  output*  but  due  to  tine  liaitations  was  unable  to  ooaplete  it. 

A  push-pull  transmitter  is  more  desirable  for  an  experiment  of  this 
nature*  but  it  is  not  absolutely  necessary  if  one  is  not  available. 

Probe  QlCSttil 

Schematically  the  circuitry  for  plotting  of  probe  curves  is  as 
shown  in  Fig.  12  on  the  following  page.  Bj*  the  voltage  dropping 
resistor*  wust  ho  email  in  oonparison  to  the  apparent  Impedance  of 
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the  plasma  Rp  so  that  the  sweep  potential  will  be  the  potential  that 
appears  betveen  probe  and  plasau,  i.e. , 


*s  *  1  X  **p  +  ty  *  1  x  *p  (8*) 

This  is  necessary  so  that  for  every  valve  of  probe  current  the  cor¬ 
responding  value  of  probe  voltage  will  be  known  accurately.  Values 
of  5,  10,  and  500  okas  were 

The  voltage  across  R^ 
is  the  input  to  the  I  axis 
of  an  X-I  plotter  and  the 
voltage  of  the  variable 
sweep  supply  is  the  input 
to  the  X  axis  of  the  plot¬ 
ter.  As  the  power  supply 
is  swept  manually  through 
various  potentials  of  the 
probe  with  respect  to  the 
reference  electrode,  the 
plotter  will  automatically  plot  the  voltage  across  R^  (which  is  pro¬ 
portional  to  probe  current)  versus  the  probe  potential. 

In  this  investigation  an  instrument  designed  by  the  Uectronic 
Technology  Laboratory  was  used  in  taking  probe  curves.  It  not  only 
plotted  probe  curves  using  an  X-I  plotter,  but  also  had  the  capability 
of  plotting  with  the  same  X-I  plotter  the  second  derivatives  of  the 


used  for  Rj  to  meet  this  requirement. 
Y 


A 

I  Rasa  re  *«e 


Fig.  12 

Circuit  for  Automatic  Probe  Plotting 
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probe  curves;  hence,  it  will  be  referred  to  as  the  Seooud  Derivative 
Apparatus. 

To  take  second  derivatives  the  Second  Derivative  Apparatus  super- 
iaposes  a  very  snail  pure  1-kc  sinusoidal  signal  on  the  DC  probe  volt¬ 
age  which  is  swept  manually.  In  this  investigation  a  40  nv  signal  was 
used.  The  second  harmonic  in  the  probe  current  caused  by  the  non¬ 
linearity  of  the  probe  curve  is  filtered,  amplified,  rectified  by  a 
phase  true  detector,  and  plotted  by  the  X-Y  plotter.  The  second 
harmonic  is  proportional  to  the  second  derivative  of  probe  current 
with  respect  to  probe  voltage.  A  more  detailed  description  of  the 
Second  Derivative  Apparatus  may  be  found  in  an  article  published  in 
the  March,  1963,  issue  of  the  Review  of  <tM ent^la  JfeAlBIMnlA  (Ref  9). 

Special  effort  was  taken  in  grounding  and  shielding  the  test 
apparatus.  The  transmitter  chassis,  the  aluminum  box  containing  the 
T-aplice,  the  reference  electrode,  the  electrical  center  of  the 
exciting  coil,  and  the  outer  conductors  of  the  coaxial  cables  were 
all  tied  into  a  common  ground.  The  radiating  sources  of  R?  power  such 
as  the  transmitter,  the  power  leads,  and  the  T-splioe  were  shielded  so 
that  the  KF  power  fed  into  the  plasma  would  be  only  that  fed  in  by  the 
exciting  coil  and  hence  would  be  introduced  aymatr ic ally .  The  Second 
Derivative  Apparatus  was  shielded  by  its  chassis,  and  its  chassis  was 
grounded. 
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Lwr*mt 

The  layout  of  the  equipment  as  used  in  this  experiment  it  shown 
in  Fig.  13  on  the  following  page.  The  large  instrument  in  the  left 
foreground  is  the  Second  Derivative  Apparatus.  Centered  in  the  pioture 
is  the  glass  tube  which  contained  the  plaasa.  The  probe  shaft  is  seen 
protruding  fron  the  top.  Surrounding  the  tube  is  the  one-turn,  oopper 
coil  used  to  excite  the  plaama,  and  to  its  right  is  the  alumlmmi  box 
which  contained  the  T-splice  used  in  changing  the  tranaadtterb  single- 
ended  feed  to  push-pull  feed.  The  instrument  to  its  rear  with  the  large 
dial  in  the  center  and  the  smaller  dial  in  the  upper  right-hand  corner 
is  the  Viking  Ranger  RF  tranaaitter.  The  large  hexagonal  object  in 
the  center  of  the  upper  part  of  the  pioture  is  the  oven  used  to  bake 
the  tube.  Figure  14  on  page  47  is  a  close-up  of  the  tube  and  excita¬ 
tion  coil.  The  coaxial  cable  fastened  to  the  probe  and  reference 
electrode  by  the  alligator  clips  leads  tq  the  Second  Derivative 
Apparatus.  The  glass  bottle  containing  the  filling  gas  nay  be  seen 
in  the  right  rear.  Valves  marked  2  and  3  in  Fig.  10  may  be  seen  in 
the  right  foreground  standing  upright  on  the  thin  Kovar  tubing  coming 
from  the  bottle  containing  the  filling  gas. 

Before  taking  probe  curves,  two  tasks  remained  to  be  accomplished. 
First,  the  exciting  coil  had  to  be  centered  around  the  plasma  tube  so 
that  the  probe  would  be  in  the  field-free  region  and  so  that  no  RF 
potential  would  exist  between  probe  and  plasma.  Secondly,  the  probe 
had  to  be  cleaned  by  sputtering  before  taking  probe  curves. 


45 


OHty'Ffagra/63-3 


fc’ig.  13 

equipment  Layout 
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Pig.  14 

Plasma  Tube  and  Exciting  Coil 
47 
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fl£  Jfeil 

The  centering  of  the  ooil  was  verified  by  the  oiroutt  in  Vlg.  15. 
The  variable  supply 
consisted  of  batteries 
in  series  to  give  poten¬ 
tials  batwaan  tha  probe 
and  rafaranoa  alaetroda 
of  0,  11,  12,  14.5,  22, 
and  29  volt a.  Thus 
points  along  tha  proba 
curve  from  tha  ion 
region  to  the  eleotron 
accelerating  region 
could  be  examined.  A 
floating  DC  sweep  supply 
was  not  used  because  it  was  fait  that  any  RF  might  be  bled  off  to 
ground  through  the  DC  supply  before  reaching  the  10  ohm  resistor.  A 
0.004  |if  RF  by-pass  capacitor  was  placed  across  the  terminals  of  the 
battery  to  provide  a  low  impedance  path  for  the  RF.  At  the  frequency 
of  operation  the  impedance  of  the  capacitor  is  less  than  2  ohms.  The 
batteries,  the  by-pass  capacitor,  and  the  10  ohm  resistor  ware  placed 
in  an  aluminum  box,  as  indicated  by  the  dotted  line  in  Fig.  15,  to 
prevent  RF  pickup.  All  leads  into  and  out  of  the  box  ware  shielded. 

The  voltage  drop  across  the  10  ohm  resistor  was  fed  into  the  Y 
axis  of  a  Tektronix  Type  545  oscilloscope.  The  X  axis  was  swept 


r - - - 


Circuit  for  Bateotian  of  V  on  Proba 
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internally.  The  coil  was  first  csntarad  about  tha  probe  visually.  At 
potentials  of  0,  11,  and  12  volts,  no  RF  appeared  on  the  scope  with 
the  vertical  axis  set  at  its  maximum  sensitivity,  i.e.,  1  av  per  inoh. 
At  14.5  volts  a  small  but  detectable  amount  was  observed,  i.e.,  about 
0.25  mv  peak-to-peak  voltage.  When  the  coil  was  adjusted  slightly  the 
RF  was  eliminated  for  all  bias  voltages  up  to  29  volts. 

mining  of  Probe 

The  probe  was  cleaned  before  data  taking  by  sputtering.  The 
probe  was  made  very  strongly  negative  with  respect  to  the  plasma. 

The  positive  ions  thus  severely  bombarded  the  probe,  causing  the  sur¬ 
face  of  the  probe  to  be  eroded  away.  The  probe  was  sputtered  with  a 
200-volt  negative  potential  for  30  seconds  on  two  separate  occasions. 
In  addition,  in  order  to  supply  additional  electrons  to  light  the 
discharge,  a  Tesla  coil  was  placed  in  contact  with  the  probe  shaft. 
This  caused  the  probe  to  become  strongly  negative  with  respect  to  the 
plasma  so  that  each  time  the  plasma  was  ignited  the  probe  was  effec¬ 
tively  sputtered. 
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vi.  Halit  and  Qtaototitaa 


at  Ha  Oft  Mad— 

The  plaana  had  two  stead j  state  aodu  of  operation.  Id  what  will 
ha  eallad  tha  waak  plwa  tha  percentage  of  ionisation  waa  ao  low  that 
a  500-obn  voltage  dropping  rasiator  waa  uaad  in  ardar  to  taka  a  proba 
curve.  Tha  value  of  tha  proba  currant  drawn  at  pi  sans  potantial  for 
tha  proba  6.3  on  fron  tha  rafaranoa  alactroda  waa  11  nioroanps.  At 
thia  position  plaana  potantial  ooourrad  at  about  22.5  volts  abova 
ground  potantial ,  Fig.  16  as  shown  on  tha  following  paga.  Using 
*q  (59) 

22.5  *  11  x  ICf6  X  (1^  ♦  500)  (60) 

R_  =  2.2  x  106  -  500  x  2.0  x  106  ohns  (61) 

P 

Thus,  tha  voltaga  dropping  rasiator  is  snail  cohered  to  tha  iapadanoa 
offarad  tha  circuit  by  tha  plaena,  and,  therefore,  tha  awaap  potantial 
(in  this  oaaa  22.5  volts)  is  tbs  potantial  that  appears  batwaan  tha 
proba  and  plaana.  Unfortunately,  tha  currant  drawn  for  tha  waak  plaana 
was  too  snail  a  value  for  tha  Saeond  Derivative  Apparatus  and  no  saoond 
derivative  curves  could  be  taken  in  this  node  of  operation. 
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v  - *- 

(  with  rtftnw 


rig.  16 

Probo  Cnrvo  of  Dla  Mod* 
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The  Moond  sod*  of  operation  of  tba  plasaa  was  considerably 
brighter  in  appear anoa  indicating  a  higher  degree  of  ionisation  and 
excitation  than  in  the  weak  node.  The  denser  node  was  achieved  by 
putting  a  Tesla  coil  in  contact  with  the  probe  shaft  while  the  plasaa 
was  in  the  weak  node.  This  produced  an  additional  source  of  electrons 
which  in  turn  increased  ionisation. 

Plasaa  potential  at  6.3  oa  froa  the  referenee  electrode  for  this 
■ode  occurred  at  about  14  volts  above  ground  potential  (Pig.  17  on  the 
following  page).  The  current  drawn  at  plasaa  potential  was  1 .4  ailli- 
aaps.  Since  the  current  was  relatively  large,  a  5-oha  voltage  dropping 
resistor  was  used  in  plotting  the  probe  curve  and  the  Second  Derivative 
Apparatus  was  able  to  take  the  second  derivative  of  the  probe  curve. 

The  second  derivative  curve  is  shown  in  Pig.  18  on  page  54. 

A  rough  idea  of  how  the  energy  distribution  differs  in  the  two 
■odes  of  operation  nay  be  obtained  by  exaaining  the  two  probe  curves 
without  recourse  to  their  second  derivatives.  The  fact  that  plasaa 
potential  is  at  a  higher  potential  with  respect  to  the  reference  elec¬ 
trode  in  the  weak  plasaa  than  in  the  dense  plasaa  suggests  that  higher 
energy  electrons  are  collected  by  the  probe  in  the  weak  plasaa  than 
in  the  dense  plasaa. 

Pron  this  fact  it  is  logical  to  asaune  that  in  the  vicinity  of 
the  probe  the  electric  field  is  stronger  in  the  weak  plasaa  than  in 
the  dense  plasaa  since  the  local  value  of  the  electron  energy  is 
fixed  by  the  local  value  of  S  (Ref  3:2784).  This  would  suggest  then 
that  as  the  plasaa  becaaes  more  ionised  it  has  a  tendency  to  shield 
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(  w fmA  H  tyw 

Tig*  17 

Proba  Curve  of  Dons*  Plasm 
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rig.  is 

Saeoad  DariYatlva  Carr*  for  Probo  at 
6.3  ea  froa  Jtafaraoea  Klaotroda 
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oat  the  W  field.  This  sffsot  know  as  skin  offset  has  bass  investigated 
fey  lekert  (Rsf  3i2788).  is  tbs  conductivity  of  the  pi  a— a  increases, 
tbs  radial  distribution  of  alaotrio  field  Intensity  changes  (qualita¬ 
tively)  trm  a  plot  indicated  fey  rig.  19a  to  a  plot  indicated  fey  fig.  19b. 


Fig.  19 

Distribution  of  Field  in  Fill 
_ (Fro*  Ref  3«2788) 


Probe  curves  and  their  seeond  derivatives  for  the  dense  plasma 
were  taken  for  the  probe  at  2.5,  3.8,  5.1,  6.3,  and  8.9  on  flros  the 
reference  electrode.  The  probe  did  not  nova  freely  in  the  tube;  hence 
probing  along  the  axis  of  the  tube  at  equally  spaced  Increments  was 
not  possible.  The  second  derivative  for  the  probe  2.5  on  from  the 
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reference  electrode  is  shown  in  Fig.  20a  on  the  following  page.  The 
resultant  energy  distribution  curve  and  a  Maxwellian  having  the  saae 
average  energy  are  shown  in  Fig.  20b.  Planes  potential  occurs  where 
the  second  derivative  curve  (Fig.  20a)  crosses  the  X  axis.  The  probe 
voltage  with  respect  to  the  plasma  is  aero  at  this  point.  To  the 
left  of  this  point  is  the  electron  retarding  region,  the  region  of 
interest.  The  ordinates  for  the  energy  distribution  curves  ware 
obtained  by  Multiplying  each  ordinate  of  the  second  derivative  curve 
in  the  electron  retarding  region  taken  at  0.1  ev  Intervals  tines  the 
square  root  of  the  corresponding  abscissas.  Thus  Eq  (46)  was  applied 
throughout  the  electron  retarding  region  to  obtain  the  energy  dis¬ 
tribution  spectra.  The  constant  K*  was  ignored  in  the  calculations 
because  at  this  stage  of  the  calculations,  n  is  not  known.  Thus 
the  ordinates  of  F(E)  nay  be  found  but  nay  not  be  normalized 
using  Eq  (46). 

Since  these  calculations  are  rather  tedious  when  done  by  hand 
for  several  probe  positions  in  the  tube,  a  computer  program,  which 
nay  be  found  in  Appendix  A,  was  written  to  accomplish  this.  This 
program  automatically  adjusts  the  RF  distribution  ordinates  so  that 
the  data  may  be  used  to  plot  a  normalized  distribution  curve.  The 
program  also  confutes  the  average  energy  of  the  electrons,  their 
average  velocity,  the  denaity  of  the  plasma  and  the  ordinates  of  a 
Maxwellian  Distribution  having  the  same  energy  as  the  average  energy 
of  the  RF  distribution. 
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(b) 

Inert j  Distribution  Curve 

flf.  20 

Second  Derivative  Curve  and  Energy  Distribution  Curve 
for  Probe  at  2.5  on  fron  Referenoe  Electrode 
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Plots  of  the  electron  energy  distribution  for  the  probe  st  3.8, 

3.1,  6.3,  and  8.9  cm  ftcm  the  reference  electrode  nay  be  found  in 
Fig.  21  (pegs  39)  and  Fig.  22  (page  60).  Plotted  with  each  RF  dis¬ 
tribution  curve  is  a  Maxwellian  curve  having  the  same  average  energy. 
Semi-log  plots  of  the  second  derivative  curves  nay  be  found  in 
Appendix  B.  Since  the  semi-log  plot  of  the  second  derivative  of  a 
Maxwellian  energy  distribution  function  is  linear,  these  curves  pro¬ 
vide  an  additional  means  of  examining  the  deviations  of  the  RF  elec¬ 
tron  energy  distribution  from  a  Maxwellian. 

In  examining  Figs.  20b,  21,  and  22  several  facts  beoome  apparent. 
First,  the  energy  distribution  curves,  although  varying  in  average 
energy,  have  generally  the  same  shape  for  the  probe  at  all  positions 
except  at  2.3  cm.  The  proximity  of  the  probe  to  the  reference  elec¬ 
trode  at  2.3  cm  probably  accounts  for  the  distortion  noted  at  this 
position.  Secondly,  the  curves  deviate  markedly  from  Maxwellian  curves. 
This  is  to  be  expected  since  a  Maxwellian  is  only  anticipated  when  the 
particles  are  in  thermal  equilibrium  and  there  is  no  external  field 
(Ref  6:24).  The  low  energy  electrons  found  in  Maxwellian  distributions 
are  missing  from  the  RF  distribution.  This  is  to  be  expected  because 
the  RF  field  penetrates  the  plasma  and  removes  the  low  energy  electrons 
by  giving  them  additional  energy.  The  high  energy  electrons  are  missing 
because  in  the  vicinity  of  the  probe  the  accelerating  electric  field 
must  be  low  (Fig.  19)  and  also  because  the  higher  energy  electrons  lose 
their  energy  by  inelastic  collisions.  There  are  two  metastable  states 
of  argon  occurring  at  11.3  and  11.7  ev  (Ref  16:39).  Thus  electrons 
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(*) 

Probe  at  3.8  cm  froa  Reference  Electrode 


C V - ► 

(b) 

Probe  at  5.1  on  froa  Reference  Electrode 


Pig.  21 

Energy  Distribution  in  RP  Plasaa 
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Proto  at  6.3  m  from  Rafsraooa  Blaetroda 


(b) 

Proto  at  8.9  oa  froa  Rafaranca  Ilaotroda 


Pig.  22 

Enargy  Distribution  in  RF  Plana 


(2lVPhys/63-3 

with  those  energies  or  higher  can  excite  argon  atoms  and  hence  will 
lose  some  or  all  of  their  energy.  The  first  ionization  potential  of 
argon  occurs  at  15.68  ev  (Ref  5:2549).  As  the  energy  distribution 
curves  in  Fig.  11  indicate,  no  electrons  collected  by  the  probe  at 
any  of  the  points  vhere  probe  measurements  were  made  were  energetic 
enough  to  cause  ionization.  At  least  they  were  not  present  in 
detectable  quantities. 

The  fact  that  the  discharge  was  self-sustaining  may  be  explained 

by  two  possibilities.  First,  it  has  been  mentioned  that  argon  has 

two  metastable  states.  Metastable  states  have  the  characteristic  of 

being  cooperatively  long-lived.  Their  decay  time  is  on  the  order  of 

—8 

milliseconds  versus  a  decay  time  of  10“  seconds  for  a  non-metastable 
state  (Ref  14*15).  Thus  it  is  possible  for  ionisation  in  the  tube  to 
occur  in  steps,  i.e.,  an  argon  atom  may  be  excited  by  an  electron  and 
then  subsequently  ionized  by  another  electron. 

Another  possibility  is  that  ionization  takes  place  at  the  outer 
extremities  of  the  tube.  This  is  quite  logical  because  as  mentioned 
previously  the  local  value  of  the  electron  energy  is  fixed  by  the 
local  value  of  E.  Since  the  field  intensity  decreases  with  the 
radius  as  Fig.  19  indicates,  the  electron  energies  would  be  expected 
to  have  lower  values  at  the  center  of  the  tube  than  at  its  outer 
extremities.  Diffusion  into  the  center  of  the  tube  by  the  fast  elec¬ 
trons  is  hindered  by  the  electric  field.  The  electrons  will  have  a 
tendency  to  follow  the  £  field  lines  which  are  radial  about  the 
center  of  the  tube.  Electrons  penetrating  into  the  center  would  have 
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to  bo  scattered  inwardly  and  this  inward  aeattor  would  bo  hiadorod  by 
inolastle  oeattorlng.  Tho  fiold  should  have  a  significant  offoet  on 
tho  notion  of  tho  olootrons  because,  as  will  bo  seen,  tho  ratio  of 
collision  frequency  to  fiold  frequency  was  about  3  to  1 . 

Table  I  shows  the  variation  of  tho  average  energy  of  tho  oloo¬ 
trons  and  also  tho  variation  of  plaaaa  density  for  various  points 
along  tho  axis  of  tho  tube.  Many  additional  ■eaonrsaonts  along  the 
axis  of  tho  tuba  need  to  bo  node  before  an  oxpltatioti  of  this  varia¬ 
tion  is  attempted. 


Table  I 

Average  Xnergy  of  Xleetrons  and  Plasm  Density 
at  Various  Probe  Positions  Along  Tube  Axis 


Position  Average  Xnergy  Plaaaa  Donoit 


lev) 

2.5 

4.67 

3.08  x  1010 

3.8 

5.13 

3.13  x  1010 

5.1 

4.84 

3.48  x  1010 

6.3 

4.84 

3.28  x  1010 

8.9 

4.84 

3.0  x  1010 
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Effect  fi£  Maenatle  rtald 

It  is  obvious  that  tha  R7  slactric  field  prevents  the  electrons 
fro*  achieving  a  Maxwellian  distribution  by  continually  supplying 
energy  to  then.  What  effect,  if  any,  a  nagnetic  field  has  on  the 
electrons  was  investigated  by  supplying  an  external  DC  nagnetic  field 
to  the  plasma. 

From  a  classical  physics  standpoint,  a  nagnetic  field  does  not 
change  the  magnitude  of  the  velocity  of  the  charged  particle  but  only 
its  direction.  Thus  the  kinetic  energy  of  the  particle  should  remain 
the  sane.  Probe  curves  were  taken  with  and  without  an  external  mag¬ 
netic  field  present,  and  they  indicated  that  the  density  of  the 
plasaa  was  increased  by  a  factor  of  1 .6  with  the  nagnetic  field  preaent. 
Second  derivative  curves  were  also  taken,  and  the  second  derivative 
curve  with  the  nagnetic  field  present  matched  perfectly  the  curve 
without  the  magnetic  field,  Indicating  that  this  external  nagnetic 
field  had  no  effect  on  the  energy  distribution  of  the  electrons.  Thus 
it  is  reasonable  to  assume  that  the  alternating  nagnetic  field  present 
in  the  plasma  does  not  affect  the  energy  distribution  of  the  elec¬ 
trons  either. 

Deriving  a  distribution  function  that  would  predict  the  energy 
distribution  of  the  electrons  in  the  tube  is  an  exceedingly  difficult 
task  and  was  not  undertaken  by  this  author.  It  is  complicated  by  the 
fact  that  the  distribution  function  will  be  a  function  of  the  electric 
field  intensity  and  this  intensity  varies  with  the  radius  of  tha  tube. 
The  problem  is  even  further  complicated  by  the  fact  that  in  this 
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Investigation  the  collision  frequency  as  will  be  seen  in  the  following 
paragraph  was  comparable  to  the  exciting  frequency.  Thus,  the  task  of 
constructing  a  model  which  describes  the  events  occurring  in  the  tube 
is  vary  difficult. 

P-Vlrml,-Mn«  £  GPIUpIQB 

The  calculation  of  the  electron  collision  frequency  in  the  danse 
•ode  was  done  by  first  determining  the  electron  mean  free  path  for 
elastic  collision  by  use  of  the  following  formula  (Kef  2:3)  : 


X.  » 


(62) 


where  X  is  mean  free  path, 

PQ  is  a  measure  of  collision  probability,  and 

PQ  is  the  reduced  pressure  in  m  Mercury,  i.e.,  PQ  =  22Z2. 
u  u  T 

Using  4.8  ev  as  the  average  energy  of  the  electron  PQ  =  32  (Kef  2:7) . 

With  T  =  293°K  Eq  (62)  becomes 

X  = - 2222s - =  1 .58  cm 

.2  x  273  x  32 

Collision  frequency  is  related  to  X  by  the  following  equation 

v  =  I  (63) 
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Using  X  *  1 .98  and  1 .16  *  107  o^seeond  as  ths  value  of  tba  average 
elactron  velocity  which  was  calculated  by  tba  computer 

V  =  lalj  S-lfl-  =  73.4  x  106  collisions/sec  (64) 

I  *Do 


Oala&aSiAm  fi£  Escort  s£  T«»t— ««■» 

Tba  percentage  of  ionisation  of  tba  gas  in  tba  bright  node  was 

Q 

also  calculated  by  first  finding  tba  nunber  of  gas  atoms/cn  in  tba 
tube  before  the  generation  of  tbe  plana  and  then  comparing  this  auabar 
to  tba  plana  density  which  was  calculated  by  the  oonputer.  Letting 
the  volune  of  the  gas  in  the  tube  be  V,  at  0.2  m  Mercury  at  293°K, 
tbe  volune  of  the  gas  at  8TP  is  given  by 


.2  x  V.  x  273 

Vn  = _ .1, _  =  .0002497. 

0  293  x  760  1 


(65) 


23 

Since  there  are  6.025  x  10  atons/nole  and  a  mole  of  argon 

3 

occupies  22,400  on  ,  the  nunber  of  argon  atoms  in  the  tube  is 


6.Q25  g 


22,400 


-lfl^x 


,000245V.,  =  6.59V1 


x  1015  atone 


(66) 


6.59V1  x  10 


,19 


=  6.59  x  1015  atoms/cm' 


(67) 


The  average  of  the  plasma  density  in  Table  I  is  3.31  x  10^.  Therefore 
for  the  dense  node  along  the  tube  axis 


%  ionisation  =  3*31  x,  .10^  x  100  =  .00050% 
6.59  x  1015 


(68) 
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Thesis  at  team 

Th*  principal  sources  of  calculable  error*  in  the  experiment  may 
be  attributed  to  reeding  of  the  ordinates  of  the  second  derivative 
curves.  The  values  of  the  average  energy  are  only  good  to  two  signifi¬ 
cant  figures.  Different  data  frcsi  th*  sane  second  derivative  curve 
were  fed  into  the  coaputer  on  two  different  occasions  and  the  dif¬ 
ference  in  average  energy  was  0.03  ev.  The  same  thing  was  tried  on 
another  curve  giving  a  difference  in  average  energy  of  0.08  ev. 

Another  source  of  error  is  in  the  probe  itself.  The  probe  does 
disturb  the  plasma  to  a  certain  extent.  Studies  made  in  the  Electronic 
Technology  Laboratory  indicate  that  the  distortion  is  most  apparent  at 
potentials  near  plasma  potential.  This  would  introduce  errors  in  the 
low  end  of  the  distribution  spectrum.  The  amount  of  distortion  intro¬ 
duced  by  the  probe  has  not  yet  been  fully  investigated. 

In  deriving  Eq  (46)  it  was  assumed  that  the  current  collected  by 
the  probe  was  purely  an  electron  current.  In  actuality  in  the  electron- 
retarding  region  the  current  is  composed  of  both  electrons  and  ions. 

The  ion  current  is  approximately  linear  in  the  ion  saturation  region 
as  Fig.  17  indicates  and  may  be  linearly  extrapolated  to  zero  for  less 
negative  voltages  (Ref  7j594) .  If  it  were  completely  linear  in  the 
extrapolated  region,  then  the  ion  current  would  have  no  effect  on  the 
second  derivative  curve  as  the  second  derivative  of  a  linear  equation 
is  everywhere  zero.  If  it  is  not  linear,  the  effect  will  still  be 
small  because  the  ion  current  must  monotonically  go  to  zero  at  a 
potential  slightly  more  positive  than  plasma  potential  and  in 
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wining  Fig.  17  this  doss  not  permit  ths  ion  current  to  be  very  con¬ 
cave  downward.  Thus  any  errors  introduced  by  the  ion  current  in  the 
second  derivative  curves  end  in  the  resultant  energy  distribution 
curves  should  be  extremely  snail. 

CMBtoalgM  and  Binn— inlitlnnn 

The  objective  of  determining  the  feasibility  of  wiring  probe 
energy  distribution  measurements  in  an  RF  plaema  was  accomplished 
and  energy  distribution  measurements  were  made  along  the  axis  of  the 
tube.  The  results  indicated  that  the  assumption  (which  is  often 
made  for  convenience)  that  the  distribution  is  Maxwellian  is  not  a 
good  one. 

There  are  many  additional  experiments  which  should  be  performed 
on  the  tube  in  order  to  understand  more  fully  the  electron  energy 
distribution  of  an  RF  plasma.  Measurements  should  be  made  at  dif¬ 
ferent  tube  pressures  to  see  what  effect  changing  the  H/p  ratio  has 
on  the  energy  distribution.  The  exciting  frequency  should  also  be 
varied  to  see  what  effect  it  has  on  the  plasma  energy  distribution, 
if  any.  This  latter  experiment  requires  a  push-pull  transmitter 
since  variable  frequency  measurements  with  the  experimental  set-up 
used  for  this  report  are  Impractical. 

in  experiment  should  be  devised  to  see  how  isotropic  the  electrons 
are  along  the  axis  of  the  tube.  This  could  be  accceaplished  by  changing 
the  orientation  of  a  plane  probe  in  the  tube  and  examining  the  resulting 
probe  curves.  This  experiment  is  Important  because  spherical  probe 
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theory  ie  beeed  on  ieotropio  distribution  of  electrons.  It  the  field 
la  weak  near  the  ads  of  tha  tuba,  the  distribution  nay  very  well  be 
quite  isotropic,  but  this  should  be  investigated. 

If  the  tube  is  constructed  with  a  snail  Movable  loop  that  oac  be 
extended  fron  the  tube  wall  to  the  osntsr  of  tha  tube,  than  tha  mag¬ 
netic  field  nay  be  probed.  If  the  magnetic  field  is  saoro  in  the 
Interior  of  the  plasma,  then  the  aeaoeiated  electric  field  must  be 
sero  and  the  plaaaa  is  thus  effectively  acting  as  a  shield  to  the 
electric  field. 

Finally,  since  this  experiment  is  so  dependent  on  symstij  it 
is  recosnanded  that  in  future  experiments  a  reference  electrode 
identical  to  the  one  used  in  this  experiment  be  placed  in  the  bottom 
of  the  tube  to  further  sysnetrise  the  conditions  in  the  plasma. 
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Appendix  A 

Caapit«r  £cSgCM  far  Data  Redaction* 

The  pro  gran  used  in  reducing  the  date  usee  lq  (46)  to  obtain  the 

unnormallsed  energy  distribution  ordinates.  The  Integral  P5®**  p(g)<jE 

% 

is  evaluated  and  nonullsed  to  1  by  Multiplying  it  by  a  normalising 
factor.  Simpson's  Rule  is  used  in  performing  the  integration.  Then 
each  energy  distribution  ordinate  calculated  by  Iq  (46)  is  multiplied 
by  the  seme  normalising  factor  mentioned  above  and  thus  the  distribu¬ 
tion  curves  plotted  on  pages  57,  59,  and  60  are  normalised. 

This  program  uses  the  AFU  Fortran  System.  The  AP1T  Processer, 
Loader  and  Subroutine  decks  nay  be  obtained  from  the  IBM  Library. 

For  input  data  an  even  number  of  equal  increments  mist  be  taken  along 
the  horizontal  axis  of  the  second  derivative  curves  since  as  mentioned 
above  Simpson's  Rule  is  used  in  the  normalisation  of  the  energy  dis¬ 
tribution  ordinates.  The  data  must  be  read  into  the  computer  in  the 
following  order:  H,  N,  CPP,  E(I)'s,  and  SD(l)'s.  These  variables 
are  defined  in  the  Dictionary  of  Variables  shown  on  the  following 
page. 


I 


* 


* Program  for  the  IBM  1620  Digital  Computer. 


71 


GNE/Phy s/63-3 

DIMENSION  E ( 200 ) ,  SD  ( 200 ) , GE ( 200 ) 

READ, H, N , CPP 

BK-.B617E-04 

PI-3.14159 

RP-.00027 

EM-.91083E-30 

Q-. 16021  E-lB 

DO  1  1-1, N 

1  READ,E(I ) 

DO  2  1-1, N 

2  READ, SD(  I  ) 

DO  3  1-1, N 

3  GE( I )-S0( I )*SQRTf  E( I ) ) 

MM-1 

NN-N-1 

4  A-0. 

B-0. 

J-0 

5  J-J+1 
l.-2*J 
M-2*J+1 
B-B+GE (L ) 

IF(NN-l.)7,7,6 

6  A— A+GE (M) 

GO  TO  5 
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7  60  T0  (8, 1 1  ) 

8  AREE-fH*(GF.(  1  )+GE( N)+4.*B+2.*A) )/3« 
F0NE-1./AREE 

L-0 

PRINT  61 
PRINT  70 
D0  10  1-1, N 
GEf I )-GE( I )*F0NE 
PRINT  62,I,GE(I) 

L-L+1 

IF(L-3)10,9, 10 

9  L-0 
PRINT  60 

10  CONTINUE 
G0  T0  12 

1 1  AVEN- (H*(GE( 1 )+GE( N)+4.*B+2.*A) )/3. 
G0  T0  14 

12  00  13  1-1, N 

13  GE(I)-GEM)*EM) 

MM- 2 

G0  T0  4 

14  T»(AVEN*2.)/(3.*BK) 

L-0 

PRINT  63 
PRINT  70 
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00  16  1-1, N 

BB- 2 . / ( SORT  f  P I ) * ( f  BK*T )** 1 . 5 ) ) 

SO( I )»BB*SQRT(E( I ) )*EXP(-(E( « )/(BK*T))) 
PRINT  64, I , SD(I ) 

L-L+1 

IF  (1.-3)  16, 15, 16 

5  L-0 

PRINT  60 

16  CONTINUE 

00  17  1-1, N 

S0( I )-SORT( 2.*E ( I )*Q/EM) 

17  GE ( I  )«GE( I )/E( I ) 

JJ-1 

18  A-0. 

B-0, 

LI.-M-1 

D0  19  1-1, LL 

B-.5*(SD( l+1)-SD(l ))*(GEf l+1)+GE(l)) 

1 9  A-A+B 

G0  T0  (20,21  ),J«J 

20  AREV-A 
G0  T0  22 

21  AREVW-A 
AVEV-AREVW/AREV 
G0  T0  24 
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22  D0  23  1-1, M 

23  GF.f  I  )-GE(  I  )*S0(  I  ) 

JJ-2 

GO  T0  18 

PD-(CPP)/(PI*fRP**2.)*Q*AVEV)  » 

PRINT  65.AVEN 

PRINT  66,AVEV 

PRINT  67, T 

PRINT  68, PD 

60  F0RMAT(/7X) 

61  F0RMAT (/ 30X, 25HRF  DISTRIBUTION  ORDINATES) 

62  F0RMAT(2HGEl4,1H«F6.4,7X) 

63  FORMAT (/30X,33UMAXWELU AN  DISTRIBUTION  ORDINATES) 

64  F0RMAT(3HGEMI4,1H-F6.4,6X) 

65  FORMAT w7  /  20X,24HAVERAGE  ELECTRON  F.NERGY-F6.3, 1X.2HEV)  * 

66  FORMAT (///20X,26HAVERAGE  ELECTRON  VEL0C I TY-F9.0, 1 X.3HMPS) 

67  FORMAT {///20X,21 HELECTR0N  TEMPERATURE-F10. 1 , 1 X.9HDEGREES  K) 

68  F0RMAT(///2OX,15HPLASMA  DENSI TY-E14.8, IX, 17HELECTR0NS  PER  MTR) 

70  FORMAT^ / / / 7  X ) 

STOP 

END 


f 
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Appendix  B 

ta  of  Second  Derivative  Curves 
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iSiiiiliiiiiiiiliiiiiiliiiiSSBBSiiiiiiil&llSiiiiiiiiiiiili^iiliS 
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agiiglMteiE555iag3aaas§ii»agiais^iiii§igiiiiSiiSgaiiiiiiiiiik^iiii 

BaBSfMSiiSMifsssBatisss^gsiiiisasagsssssggssaasBsS-s^SBSgiiasBBSsssiag 

13SIKSin&53gl£33iiSlgif-j~ag3asC31i3£Sigfl2flNliiiaiIflI&i3ag3i£Ili33.IBS 

iSHiiilisiiSiililiilif^iiilSSiiliiiiiiiliiiiaiiiiiiiiiiiiiiiiifSSMii 
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